Optical coupling of fundamental whispering gallery modes in bi-spheres 
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What will happen if two identical microspheres, with fundamental whispering gallery modes 
excited in each of them, become optically coupled? Conventional wisdom based on coupled-mode 
arguments says that two new modes, bonding and anti-bonding, with two split frequencies would 
be formed. In this paper we demonstrate, using exact multi-sphere Mie theory, that in reality an 
attempt to couple two fundamental modes of microspheres would result in a complex multi-resonance 
optical response with the field distribution significantly deviating from predictions of coupled-mode 
type theories. 
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Introduction Optical micro-resonators attract a great 
deal of interest because of their unique properties such 
as small mode volumes, high Q-factors, etc [1]. Individ- 
ual resonators of various shapes such as microspheres 0, 
microdisks Q, microtoroidsjj], etc. have been carefully 
studied and shown to have practically achievable Q- 
factors as large as 10 11 — 10 13 . Optical coupling between 
high-Q resonators due to their evanescent fields draws at- 
tention because of unusual optical effects that can arise 
in such systems and also from the point of view of their 
potential for integrated optical circuits. Simplest of such 
configurations is a double-resonator photonic molecule 
suggested in[{| and studied both experimentally and the- 
oretically in many papers 0, 0, 

In the case of microspheres, resonances are whispering 
gallery modes (WGM) \l, m, s) characterized by angular, 
I, azimuthal, to, and radial, s quantum numbers. High 
values of Q usually corresponds to modes with I 3> 1, 
which, therefore, possess a high degree of degeneracy. 
The modes with the same orbital and radial but differ- 
ent azimuthal numbers have the same (complex-valued) 
frequency, but different spatial distributions. Of greatest 
interest are so called "fundamental modes" with \m\ = I, 
s = 1, whose field is tightly concentrated in the vicinity 
of the equatorial plane and the surface of the sphere. It is 
assumed that these modes can be selectively excited by 
coupling to a tapered fiber When microspheres are 
arranged in coupled structures the question of particular 
interest is the field distribution resulting from coupling of 
fundamental modes of individual spheres. In the frame- 
work of coupled-mode theories floL 11 | or perturbation 
approaches |l2j this distribution is often described as a 
linear combination of phase-matched modes with to = I 
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in one sphere and to = — I in the other, with the cou- 
pling strength characterize d by an overlap integral of the 
respective modal functions 12, 13}. 

However, it will be shown below that such field config- 
urations are not consistent with the symmetry of the sys- 
tem if the deviations of the spheres from the ideal shape 
is negligible. More accurately, we demonstrate, analyt- 
ically and numerically, that if one excites a fundamen- 
tal mode \L, L, 1) in one or both uncoupled spheres and 
bring them together for optical coupling, the ensuing vi- 
olation of the complete spherical symmetry of the system 
results in a field configuration containing linear combina- 
tion of all initially degenerate modes with \m\ < L. The 
optical coupling removes the degeneracy of these modes 
so that each one of them gives rise to a couple of reso- 
nance peaks at its own frequencies producing a compli- 
cated multi-peak optical response. Moreover, because of 
the coupling between modes with different I the field of 
the bi-sphere contains a significant admixture of modes 
\l, m, s) with I ^ L and s > 1. 

Fundamental modes and coordinate systems. The sys- 
tem considered in the paper consists of two identical di- 
electric spheres of radius R and refractive index n po- 
sitioned at a distance d > 2R between their centers 
(see Fig [I}. We assume that the field of the system is 
monochromatic with frequency u> and separate it into a 
combination of incident, scattered, and internal fields. 
Using a standard multi-sphere Mie theory^, [3, HH we 
present all these fields as linear combinations of single- 
sphere vector spherical harmonics (VSH). In this repre- 
sentation, the fields are characterized by expansion co- 
efficients (j oj'^, cfL for incident, scattered, and in- 
to 

tcrnal fields of TM polarization, respectively, and /„, 

^lm' *4*m f° r polarization. Here index i enumerates 
spheres and takes only two values i — 1,2. Using Maxwell 
boundary conditions and the addition theorem for the 
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VSH[H[l7| one can derive a system of equations relat- 
ing scattering coefficients a] m and b\ m to the coefficients 
of the incident field Q and rft m : 

a\,m = a ( x )\ N) I 0,m + E E [ a l',m' A l,m ( s > r i - r 

)]} 



A z, X 



u l',m' D l,m \ 



(1) 



where r 

„(M) 



is the position vector of i-th sphere, (ty and 
a; ' are single sphere Mie scattering coefficients for TM 
and TE polarizations respectively. These coefficients 
have poles at the specific values of the dimensionless fre- 
quency parameter x, defined as x = nRuj/c, which de- 
termine frequency and the spectral width of the single- 
sphere WGM resonances. Explicit expressions for the 
scattering coefficients as well as for translational coeffi- 
cients A l '™ {vj — Yi) and B\ (rj — n), which describe 
optical coupling between the spheres can be found, for in- 
stance in Ref. [1. [13. flit . 

Our goal is to study effects of optical coupling on 
fundamental modes of single spheres. Classification of 
WGMs as fundamental, however, is linked to a particu- 
lar coordinate system: modes with |m| = I are associated 
with the plane perpendicular to polar axis z. In order to 
achieve strong coupling between fundamental modes in a 
bi-sphere one needs to choose a polar axes such that the 
respective equatorial plane would contain centers of both 
spheres. This coordinate system is labeled by lower case 
letters in FigfT] However, this coordinate system is not 
consistent with the symmetry of the bi-sphere, so that its 
modes defined using associated spherical coordinates can- 
not be characterized by azimuthal number m. Formally, 
it follows from the presence of non-diagonal in to elements 
in the translation coefficients. Consequently, the config- 
uration of the fields in this system cannot be presented 
as a combination of the modes with \m\ = I contrary to 
the assumption of the coupled-mode theories. In a coor- 
dinate system with polar axis along the line of symmetry 
of the system (XY Z system in FigfT]) the translation co- 
efficients are diagonal in m and the normal modes of the 
coupled system can again be classified according to the 
azimuthal number 0, However, the field distribution 
corresponding to the fundamental mode of the xyz coor- 
dinate system cannot be characterized by a VSH with a 
single m in the XYZ system. In order to describe the 
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FIG. 1: Configuration of a two-sphere system and possible 
coordinate systems. Arrowed lines show schematically two 
single-sphere fundamental modes 



same field distribution in the new coordinates one has to 
use transformation properties of VSH 15] and express a 
single \L,L, 1) mode of xyz system as a linear combina- 
tion of all \L, to, 1) modes of the XYZ system. In order 
to reproduce such field configuration in a single sphere 
we choose the incident field to be TE polarized with ex- 
pansion coefficients of the form 



5lL&ilRLm, Rhr, 



2 L 



{2L)\ 



(L + m)\(L — to)! 

(3) 

which is dictated by the transformation properties of 
VSH. 

Fundamental modes in the bi-sphere Since transla- 
tional coefficients in the new coordinate system become 
diagonal in to, expansion coefficients in Eq.(UJ and ^ 
with different values of the azimuthal number become 
independent and can be solved separately. In what fol- 
lows we will also neglect the cross-polarization terms, 
described by coefficients B\™, which are usually much 
smaller than Then we need only to be concerned 

with coefficients bi m and Eq.([2]). In order to understand 
qualitatively the effect of the coupling on the field distri- 
bution in the bi-sphere we will first solve this equation in 
so called single-mode approximation neglecting interac- 
tion between modes with different ^-numbers (see details 
in Ref. [H, [l8| ) . Expression for the scattered field in this 
approximation can be presented in the following form 



2 ^ 



Rhr, 



M Lm (r-n) + M £m (r-r 2 ) 
"^-^(ra-n) 
M iOT (r - n) - M Lm (r - r 2 ) 



(4) 



which gives a clear physical picture of the phenomenon 
under consideration. This expression describes an optical 
response with resonances at two sets of frequencies: one is 
given by the zeroes of a~^ 1 — Aj^™ anc ^ the other by the ze- 
roes of Uj^ + The role of coupling between spheres 
is described by translation coefficients A£m , which shifts 
resonant frequency from the single sphere values by dif- 
ferent amounts for different values of to. As a result 
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FIG. 2: Frequency dependence of the internal energy of the 
bi-sphere in the vicinity of single sphere resonance frequency 
x=21.46 corresponding to a TE mode |29, m, 1). (a) single- 
mode approximation; (b) numerical calculations including all 
modes with I < 29 



terms with different azimuthal numbers resonate at dif- 
ferent frequencies so that one ideally could expect 2(L+1) 
resonance peaks in the optical response of the system con- 
trary to the coupled-mode theory expectation of just two 
resonances. The number of observable peaks is smaller 
because frequency shifts for higher values of \m\ are small 
and resonances from respective terms may not be distin- 
guishable. 

One can notice that the numerators of terms in Eq. ([4J 
corresponding to different values of m have a form of sym- 
metric and antisymmetric combinations of respective sin- 
gle sphere modes, which is typical for the coupled mode 
approaches. In our consideration, however, this combi- 
nation appears for every m-component constituting the 
fundamental mode separately, rather than for the entire 
single sphere mode. 

Additional deviations from the coupled-mode approx- 
imation appear due to coupling between modes with dif- 
ferent values of the angular number Z, which is neglected 
in the single-mode approximation. This approximation 
can be improved by perturbation methods described in 
Ref.@, [HI, in this paper, however, we will rely on exact 
numerical treatment of Eq.([2]). For numerical analysis 
of Eq.Q with incident field in the form of Eq.Q we 
chose L = 29 and spectral interval in the vicinity of res- 
onance frequency of the mode with radial number s = 1 , 
2^29,1 = 21.46. We solve this equation using matrix in- 
version and including all coefficients with I < L. The 
limitations for the number of modes from above is made 
only to shorten the required time of calculation, but it 
does not affect the results significantly: we extended the 
number of included modes up to I — 48 for several val- 



ues of the dimensionless frequency x and observed good 
convergence of the results with only small quantitative 
changes resulting from contribution of the modes with 
I > L. Taking into account the coefficients with I < L, 
however, is significant because some modes with smaller 
values of I might have spectral overlap with the main 
mode I = L and their contribution can be resonantly 
enhanced [1, [HI]. Using found scattering coefficients we 
calculate expansion coefficients of the internal field di^ m : 
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x ji(x)[nxji{nx)]' - ji(nx)[xji{x)]' 



(5) 



where ji(x) is the spherical bessel function and [zji(z)]' 
means differentiation with respect to z. Knowing coeffi- 
cients di m we can find total energy of the field concen- 
trated inside spheres as a function of frequency, which is 
best suited to characterize optical response of our system 
in the spectral range of high-Q WGMs[8[. Fig [5] presents 
the spectra of the internal energy obtained by two pro- 
cedures: (a) calculations based on the single-mode ap- 
proximation, Eq.Q, and (b) multi-mode numerical cal- 
culations described above. One can see that while the 
single-mode model reproduces the multi-resonance op- 
tical response, it deviates strongly from exact numeri- 
cal calculations in positions and heights of the respective 
peaks. The most striking difference is significant lower- 
ing of the heights of the peaks when inter-mode coupling 
is taken into account, which indicates a large (three or- 
ders of magnitude) reduction of Q- factors of some of the 
resonances due to coupling to low-Q modes with I < L. 

In order to further elucidate the role of the inter-mode 
coupling we consider the coefficients of the internal field 
di, n for all values of < I < 29 and respective values of 
\m\ < I. In order to visualize the results we present pairs 
l,m as an one-dimensional array ordered according to 
the following rule (29, -29), (29, -28), . . . , (29, 0), (29, 1), 
. . . , (29, 29), (28, -28), ... and plot coefficients d lm ver- 
sus a number of the respective pair in this array. Fig. [3] 
presents the results of these calculations for I > 25 for 
one of the resonance frequencies, x = 21.415. One can 
see that the largest values of the coefficients correspond 
to the main angular number Z — 29, which as a function 
of to shows two large symmetric peaks for to = ±4. The 
positions of the peaks are determined by the dependence 
of the resonant frequency on the azimuthal number: the 
resonance at the chosen frequency results from the com- 
ponent of the internal field with \m\ — 4. In addition 
to main coefficients corresponding to I = 29 the plot 
reveals presence of other coefficients as well. The sec- 
ond largest coefficient corresponds to I = 25, which is 
in agreement with the fact that a frequency of the mode 
|25,m, 2) is almost in resonance with our main mode, 
|29,m, l)[8j. While individual contributions from modes 
with I < 29 appears to be small, their cumulative ef- 
fect is significant as evidenced by FigJ5] The role of the 
inter-mode coupling can be even more important if one 
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FIG. 3: Internal field coefficients for all values of I and m 
for Z > 25. Arrows indicate the regions corresponding to 
coefficients with given Z and varying m. 



attempts to excite other fundamental modes such as TE 
mode 1 39 , 39, 1). This mode interacts resonantly with TE 
modes, |34. m.2)[^. Hi| because of their large Q-factors 
and significant spectral overlap with 39,39,1). Devia- 
tions from single mode approximation in this case will 
be even more drastic. 

Conclusion In this paper we studied effects of opti- 
cal coupling on the fundamental modes excited in two 
spherical microresonators. We found that the optical re- 
sponse of this system under these excitation conditions 
is characterized by a rich spectrum with multiple reso- 
nances in contrast with just two peaks predicted by sim- 
plified coupled-mode type approaches. These multiple 
resonances arise due to the fact that strongly coupled 
fundamental modes of single spheres do not constitute a 
normal mode of a bi-sphere and must, therefore, be de- 
scribed by a linear combination of vector spherical har- 
monic with different azimuthal numbers. Optical cou- 
pling removes degeneracy between resonances with dif- 
ferent m, causing, therefore, different m-components of 
the field to resonate at different frequencies and produc- 
ing multi-peak spectrum. We also showed that coupling 
between modes with different angular numbers I also sig- 
nificantly affects spectrum and the field distribution in 
the bi-sphere. 

The question arises, however, how do our results agree 
with observations of bonding and anti-bonding orbitals 
with two split frequencies reported in many works? First 
of all, modes observed in many experiments were true 
normal modes of the bi-spheres, characterized by a well 
defined azimuthal number m @). These modes, however, 
are not coupled fundamental modes and results of this 
work do not apply to those experiments. 

In order to observe effects presented in this paper 
one needs to excite a fundamental mode in a single 
sphere concentrated in a plane that would contain cen- 
ters of both spheres, and measure optical response of the 
bi-sphere under the same excitation conditions. This 
might be possible with the use of tapered fiber exci- 



tation techniques §, one, however, should be aware of 
sensitivity of the multiple peak response to the strength 
of coupling and Q-factors of the modes. Increasing dis- 
tance between the spheres in our numerical simulations, 
we observed how multiple peaks collapse to form a double 
peak structure. It would be a mistake, however, to asso- 
ciate these two peaks with the bonding and anti-bonding 
states of the coupled-mode theory because these peaks 
arise due to smearing of original multiple peaks rather 
than due to splitting of a single sphere resonance. The 
width of those peaks would represent an inhomogeneous 
broadening caused by inability to resolve the fine struc- 
ture the spectrum rather than the radiative decay rate of 
respective modes. 

One can see, therefore, that the picture of two coupled 
fundamental modes, which is used widely in modeling of 
optical networks can be misleading, and we hope that this 
work will serve as a warning that one needs to exercise 
caution when analyzing experimental data and develop- 
ing theoretical models based on this picture. 
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